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AERONAUTICAL SYMBOLS 
1. FUNDAMENTAL AND DERIVED UNITS 



Symbol 



Metric 



Unit 



Symbol 



English 



Unit 



Symbol 



Length. 
Time___ 
Force. _ 



I 
t 

F 



meter 

second 

weight of one kilogram. 



s 

kg 



foot (or mile) 

second (or hour) 

weight of one pound. 



ft. (or mi.) 
sec. (or hr.) 
lb. 



Power. 
Speed _ 



kg/m/s_ 
/km/hr_. 
\m/s 



k. p. h. 
m. p. s. 



horsepower. 

mi./hr. 

ft. /sec. 



hp 

in. p. h. 

f. p. 8. 



m f Mass, = — 
9 



2. GENERAL SYMBOLS, ETC. 

TP, Weight, = mg mk 2 , Moment of inertia (indicate axis of the 
<7, Standard acceleration of gravity = 9. 806G5 radius of gyration, ~k, by proper sub- 

m/s 2 = 32.1740 ft./sec. 2 script). 

IF $, Area. 

8 m Wing area, etc. 

p, Density (mass per unit volume). G, Gap. 

Standard density of dry air, 0.12497 (kg-m -4 &, Span. 

s 2 ) at 15° C and 760 mm = 0.002378 (lb.- c, Chord length. 

ft." 4 sec. 2 ). b/c, Aspect ratio. 

Specific weight of "standard" air, 1.2255 /, Distance from C. G. to elevator hinge. 

kg/m 3 = 0.07651 lb. /ft. 3 /x. Coefficient of viscosity. 

3. AERODYNAMICAL SYMBOLS 



V, True air speed. 

Dynamic (or impact) pressure = 4 p V 2 

L 



1 



Z, Lift, absolute coefficient C L = 



D, Drag, absolute coefficient C D = 



qS 
D 



Cj Cross-wind force, absolute coefficient 
aS 

B, Resultant force. (Note that these cocffi- °V Center of P^^sure coefficient (ratio of 



7, Dihedral angle. 
VI 

p — * Reynolds Number, where I is a linear 
dimension. 

e. g., for a model airfoil 3 in. chord, 100 
mi./hr. normal pressure, 0° C: 255,000 
and at 15° C, 230,000; 
or for a model of 10 cm chord 40 m ; s, 
corresponding numbers are 299,000 and 
270,000. 



cients are twice as large as the old co- 
efficients L c , D c .) 
i Wi Angle of setting of wings (relative to thrust £ 
line). 

i t , Angle of stabilizer setting with reference to a, Angle of attack. 

thrurt line. €, Angle of downwash 



distance of C. P. from leading edge to 
chord length). 
Angle of stabilizer setting with reference 
to lower wing, = (i, — i w ). 
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EFFECT OF TURBULENCE IN WIND TUNNEL MEASUREMENTS 

By H. L. Dryden and A. M. Kuethc 



SUMMARY 

This investigation was carried out at the Bureau of 
Standards at the request of and with the financial assist- 
ance of the National Advisory Committee for Aero- 
nautics. The paper gives some quantitative measure- 
ments of wind tunnel turbulence and its effect on the air 
resistance of spheres and airship models, measurements 
made possible by the hot wire anemometer and associated 
apparatus developed at the Bureau of Standards. The 
apparatus in its original form was described in Technical 
Beport No. 820 and some modifications are presented in 
an appendix to the present paper. 

One important result of the present work is a curve by 
means of which measurements of the air resistance of 
spheres can be interpreted to give the turbulence quanti- 
tatively. Another is the definite proof that the discrep- 
ancies in the results on the N. P. L. standard airship 
models are due mainly to differences in the turbulence of 
the wind tunnels in which the tests were made. 

An attempt is made to interpret the observed results in 
terms of the boundary layer theory and for this purpose a 
brief account is given of the physical bases of this theory 
and of conceptions that have been obtained by analogy 
with the laws of floic in pipes. 

INTRODUCTION 

Early in the history of wind tunnel measurements it 
became apparent that there were large discrepancies in 
results obtained in different laboratories on some 
models. With improvements in technique, some of 
the discrepancies were removed or explained but there 
has always been a demand on the part of the designers 
of aircraft that wind tunnels be standardized. It was 
supposed that by a series of comparative tests, some 
correction factor could be determined by means of 
which measurements in a given wind tunnel could be 
reduced to some standard. 

In March 1920, the British Aeronautical Research 
Committee instituted a series of comparative tests to 
be conducted in as many as possible of the aerody- 
namic laboratories of the world. The purpose of the 
tests was stated in reference 1 as follows: 

"It was thought that such tests, in which the same 
models would be tested successively by all laboratories, 
would supply valuable information which had not 



previously been available. The aim of wind tunnel 
experimental work is to obtain reliable estimates of the 
forces which would be experienced by bodies moving at 
specified speeds through still air of infinite extent; but 
in practice it is necessary to hold the model stationary 
and to generate a flow of air past it and measurements 
made in this way are in some degree open to question 
in that the forces imposed upon the model may be 
affected (1) by the limited extent of the air stream in 
which they are placed and (2) by the turbulence which 
can never be entirely eliminated. The results must 
furthermore depend to some extent upon the methods 
adopted for connecting the models to the measuring 
apparatus. Different methods are adopted in differ- 
ent countries, and wind tunnels of varying size and 
design are employed; thus there is some uncertainty as 
to the extent to which a comparison can be made— 
e. g. between different airfoils tested in different coun- 
tries — and this uncertainty, it was thought, would be 
reduced if comparative figures were available from 
tests upon the same models." 

The tests are still in progress. They comprise the 
determination of lift, drag, and center of pressure for 
a standard airfoil model at various angles of attack 
and measurements of the resistance of two streamline 
models at zero angle of yaw. A report (reference 2) 
has been published on the tests of the airfoil model 
carried out in several American laboratories. The 
maximum deviations of the results from the mean 
values are of the order of 3 to 5 per cent and it is con- 
cluded that the agreement obtained in tests on airfoils 
depends almost entirely on the care used in making 
the tests. 

No report has yet been published on tests of the 
airship models but it is known (reference 3) that 
maximum deviations from the mean are of the order 
of 50 per cent and that the differences are probably 
ascribable largely to the differences in turbulence 
between the several wind tunnels. It was rather unfor- 
tunate that spheres were not included in the program 
of international tests, because spheres are also very 
sensitive to turbulence. 

Most experiments on the effect of turbulence in 
wind tunnel experiments have been qualitative in 
character and in fact in the case of the airship models 
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the observed effects have been attributed to turbulence 
on the basis of a process of elimination rather than on 
direct experiment. With the development of appa- 
ratus at the Bureau of Standards for the quantitative 
measurement of turbulence (reference 4), it became 
possible to study quanti tatively the effect of tur- 
bulence on the drag coefficient of models. This work 
was carried out at the Bureau of Standards with the 
cooperation and financial assistance of the National 
Advisory Committee for Aeronautics. 

BOUNDARY-LAYER THEORY 

The discussion of the experiments described in this 
paper will be phrased in the language of the boundary- 
layer theory of Prandtl (reference 5) and as there is no 
one article to which the reader may be referred for the 
necessary information, it is desirable to state briefly 
the elements of this theory. It is well known that in 
the greater part of the field of flow about any object 
at the Reynolds Numbers encountered in wind tunnel 
experiments or at greater Reynolds Numbers, the flow 
is approximately irrotational in character, and the 
dissipation of energy is negligible. The experimentally 
observed fact is that so long as we do not enter the 
eddying wake behind the body, the pressure on the open 
end of a tube placed parallel to the flow (which is a 
measure of the total energy per unit volume) remains 
constant throughout the field. The wake, extending 
downstream from the body, in many cases has a cross 
section equal to or slightly greater than the cross section 
of the body at the beginning and increases downstream, 
but in other cases (airship hulls or finely tapered struts) 
the wake is very small. At least over the upstream 
part of the body, the total-head tube may be brought 
exceedingly close to the surface without observing any 
change in its indication. The speed at the surface is 
known to be zero from the experiments of Stanton 
and his coworkers (reference 6), yet not far away 
from the surface, the speed is observed to be relatively 
high. These experiments indicate that the effect of 
viscosity (at least over the upstream parts of bodies) 
is confined to a very thin layer. Prandtl's introduc- 
tion of this hypothesis, namely, that the field may be 
divided into two regions, in one of which the effect 
of viscosity is negligible led to the so-called boundary- 
layer theory. 

It was found by a consideration of the order of 
magnitude of terms in the general equations of motion 
of a viscous fluid that if such a layer existed, its thick- 
ness must be of the order of magnitude of the square 
root of the product of kinematic viscosity and the 
distance from the nose divided by the square root of 
the air speed, i. e., for air, if the speed and distance 
from the nose are taken as 1 ft. /sec. and 1 foot, respec- 
tively, of the order of 0.15 inch. The equations finally 
arrived at for the steady flow of an incompressible fluid 
in the boundary layer along a 2-dimensional surface 



whose radius of curvature is large as compared with 
the thickness of the layer are as follows : 



bu . bu b 2 u 

by 

bu . bv _ 
bx~ h by~° 



lbp 
p bx 



(1) 

(2) 
(3) 



where u is the tangential component of the velocity, v 
the normal component, x the distance measured along 
the surface, y the distance measured normal to the 
surface, v the kinematic viscosity, and p the pressure. 
As boundary conditions we have u = 0, v = 0 at the 
boundary y = 0;u= U, the speed in the potential flow 
at y = 6, the outer edge of the layer. Equation (1) 
states that the tangential acceleration of a fluid particle 
is produced by the resultant of the forces due to the 
pressure and the forces due to viscosity. The term 

b 2 u 

omitted from the general equation is v g-g- Equation 
(3) is the equation of continuity. Equation (2) is 



that for the normal acceleration. The 



bv 

terms u^-> 
bx 



bv b 2 v b 2 v . . n , 

V by' P bx 2 ' v by 2 a PP earm S 111 ^ e general equations are 

neglected, being of the second order. Equation (2) 
states that the pressure does not vary across the bound- 
ary layer. It is therefore the same as the pressure in 
the potential flow outside the layer. Hence the pres- 
sure on the surface of the body is equal to the pressure 
in the potential flow at the outer edge of the layer. We 
may then compute from the observed pressure distribu- 
tion by means of Bernoulli's theorem the speed at the 
outer edge of the boundary layer and thus obtain all 
of the data necessary for a solution of the boundary 
layer equations. 

The exact solution of the equations of the boundary 
lajw (equations 1, 2, 3) meets with great difficulties, 
although by great labor a solution can be obtained in 
any numerical case. (Reference 7.) An exact solu- 
tion has been obtained by Blasius (references 8 and 9), 
by means of series developments, for the case of skin 
friction on a thin flat plate of infinite breadth. In 

this case ^ is negligible. The speed u increases 

asymptotically to its limiting value U and hence no 
exact value can be assigned to the thickness of the 

Vvx 
jY The drag 

coefficient, namely, the force divided by the product 
of velocity pressure, % p U 2 and area of the plate, comes 

out equal to 2.68-^/^* where a: is the length of the plate 

in the direction of the stream. The resistance of a 
given plate therefore varies as U L5 . 
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In the present paper we shall have little to do with 
the exact solutions, but we shall attempt to carry the 
discussion somewhat further by the approximation 
methods developed by Karman and Pohlhausen. 
(References 10 and 12.) These approximations are 
based on a new equation, an integral equation, obtained 
as a first integral of the differential equations previously 
given or derived directly from the principle of momen- 
tum. Because of the importance of this equation, 
known as the Karman integral relation, we shall repeat 
its derivation. 

The principle of momentum applied to steady 
motion states that the integral of the flux of momentum 
taken over any closed surface is equal to the resultant 
of the forces acting on the surface. 1 Let us apply this 
principle to a section of the boundary layer (Figure 1) 
of thickness 8 and width dx. The mass of fluid enter- 
ing per second through a small element dy in the 
left boundary of the section is equal to p u dy and the 
momentum per second brought in is p u 2 dy. The 
total momentum per second entering through the 

left boundary is therefore j p u 2 dy. The momentum 

per second leaving through the section at the right is 

p u l J* P u% dT/Jdx. Some momentum 

enters through the sloping upper boundary whose slope 
exceeds that of the local streamline. To obtain this 
momentum we note that the total mass per second 

entering at the left is J pu dy and the total mass per 

second leaving at the right is J p u dy + p u dyj 

dx. The mass per second entering through the sloping 

upper boundary is theref ore J pu d?/J dx and since 

this mass enters with velocity J7, the momentum per 
second entering is 

u [Lf! 

The time rate of increase of downstream momentum 
within the volume element is therefore 

This must equal the sum of the downstream com- 
ponents of all the forces acting on the surface of the 

element. These are p5 on the left boundary, 

1 — I d8 

+ ^ x (pS) dx on the right, p ^ dx on the upper sloping 

boundary and dx on the lower boundary. 

i A more precise statement is as follows: If any fixed closed surface be described 
within a steady stream of fluid, the time rate of increase, within the surface, of 
momentum in any given direction is equal to the sum of the components in that 
direction of all the forces acting on the fluid. If body forces, such as weight, are 
absent or negligible, as they are in the present instance, the only forces to be consid- 
ered are those acting at the surface, viz, normal forces due to hydrostatic pressure, 
and tangential forces due to viscous shearing. 



(/x is the viscosity). The total resultant force is 
therefore -S^dx-M^g^) dx and the principle of 
momentum slates 



d 
dx 



dp 
dx 



(4) 




I I'.ure 1.— Forces on an element of the boundary layer 

We now introduce a function q such that u=U-q*. 
We find 



d_ 



+ P l[)o <lldV l 



Hence 



^^-4?[j:H^d d ,[j:^] 



dp /dq\ 
^dx+^Vdl^o 



* q is the amount by which the longitudinal speed at any distance y from the 
solid surface is less than the speed of the free steam outside, or the retardation 
due to the proximity of the surface. 
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By Bernoulli's principle,^ (p + ^p t/ 3 )= -jf + P V ^ 

= 0, and hence the multipliers of 8 on the two sides are 

equal. We have left, dividing by p, the Karman 
integral relation, 

To use this equation for the development of ap- 
proximate solutions, any reasonable assumption is 
made as to the variation of q with y and the integral 
relation is used to give a differential equation for the 
thickness of the boundary layer. For example if we 
make the rough assumption of a linear distribution, 

«-"0 -1h H>-?- J>- ?•(!)... 

= - y and the equation for 5 is found to be 



or 



U 2 d8 5U8dU_ vU 
~ 6 dx 6 dx 8 

dU 

dx U U 



(6) 



The solution of this equation is readily seen to be 

&2 =W° f™* (7) 

and the force, F x , per unit breadth across the stream, is 
given by 



ZPd x 

-vfe/vr 

'X' 



Udx 
8 



(8) 



If we set U»dx = I, the force coefficient, C7,=rTrr- 
is found to be 



_ J\lip C* U 6 dx 



(9) 



For the case of skin friction on a flat plate, U= con- 
stant = U„. Hence /= U 0 9 x 



g2== 12vx 



U 0 



(10) 

as compared to the exact solution C F = 2.68<*U%— 

V UqX 



Still closer approximations may be made by assum- 
ing that q = a + by + cy 2 + dy* ) etc., determining the 
coefficients by suitably chosen boundary conditions, 

for example u= U at y = 8, u = 0 at y = 0, ^ = 0 at y = 8, 

2 = ~ g| at 2/ = 0 from equation (1). Enough 

boundary conditions are chosen to fix the coefficients 
in terms of 5. The integral relation then gives a 
differential equation from which 8 may be determined 
as a function of x, the distance downstream from the 
nose. The reader is referred to Pohlhausen's paper 
(reference 10) for further examples. 

SEPARATION 

When the pressure increases downstream, an inter- 
esting result is obtained by the use of four or more 

terms in the expression for q. It is found that ^ 

may vanish and the flow near the wall reverse. The 
fluid particles near the wall are dragged along by the 
friction of the neighboring particles but are retarded 
by the pressure. As the boundary layer thickens the 
retarding effect prevails and actually causes a reversed 
flow. Such reversal entails separation of the flow 
from the surface observed on cylinders, and on 
airfoils at the burble point. The boundary layer 
theory thus accounts for the breaking away of the flow 
from the surface and states that the separation is 
determined by the dissipation of energy occurring in 
the boundary layer. 

EDDYING BOUNDARY LAYERS 

The remarks made so far apply to boundary layers 
in which the flow is laminar. Experiment leads us to 
believe that the flow is more often eddying 2 as in the 
case of flow in pipes and the laws of laminar motion 
do not apply. The experimental results on flow in 
pipes are assumed to apply to the boundary layer. A 
general account of the phenomena in pipes which are 
of interest in this connection is given by L. Schiller. 
(Reference 13.) We may distinguish two values of the 
Reynolds Number (product of mean speed by the 
diameter of the pipe divided by the kinematic vis- 
cosity), namely a lower Reynolds Number, below 
which any turbulence initially present is finally damped 

2 The word turbulent is commonly used in this connection. It is desirable to 
distinguish this turbulence from the turbulence encountered in wind tunnel air 
streams. The difference is principally in order of magnitude although the tur- 
bulence in a wind-tunnel air stream is imposed by a honeycomb or other means on 
the flow from without, whereas eddying flow, as we shall term it, arises from an 
internal instability. The distinction here made will be appreciated by readers who 
have seen a demonstration of Reynolds original experiment with streams of color 
In eddying flow the stream of color is very rapidly diffused throughout the whole 
tube. The turbulence in wind-tunnel air streams corresponds to a wavering or 
fluctuation of the line of color, a turbulence of a different order of magnitude from 
that in the eddying flow and imposed from without. 

Some authors make the distinction by using" disturbance" or "initial disturb- 
ance" where we use "turbulence," and 'turbulent" where we use "eddying," 
but we feel that the use of the word turbulence to describe the departures of wind- 
tunnel air streams from uniformity and steadiness is well established. Although the 
use of "turbulence" in both cases does not in general cause contusion, we have 
thought it preferable in the interest of clarity to use different terms. 
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out and the flow becomes laminar, and an upper Rey- 
nolds Number at which laminar flow changes to 
eddying flow. The value of the low Reynolds Number 
is approximately 2,000 (Schillers value is 2,320). The 
value of the upper Reynolds Number depends on the 
amount of initial turbulence. The turbulence may 
enter the pipe in the body of the fluid or may be set up 
by the form of entrance. The highest value of the 
Reynolds Number at which laminar flow has been 
observed is 51,000 (Reference 14), although in most 
experiments values of the order of 12,000 are the highest 
observed. The results are summarized by Schiller as 
follows : 

"The result of the present work, in which the critical 
Reynolds Number is a function of the 'greatest' 
disturbance present, gives the following picture of the 
stability relations of laminar and eddying flow in 
smooth tubes. To every Reynolds Number above 
2,320 there corresponds a quite definite amount of 
disturbance which is required to produce eddying 
flow. The higher the Reynolds Number, the smaller 
is the necessary disturbance. Against still smaller 
disturbances, the laminar flow is stable. Below the 
lower critical Reynolds Number, R = 2,320, the laminar 
flow is stable against any disturbance, however great. 
There, no eddying state of flow is possible; vortices 
present will always disappear if given sufficient time." 

We shall assume the results obtained for pipes to 
apply to the boundary layer, the thickness of the 
" boundary layer" in the pipe being the radius of the 
pipe, and the speed at the outer edge of the boundary 
layer corresponding to the speed at the center of the 
pipe. Since in laminar flow in pipes, the speed at the 
center is twice the mean speed, the Reynolds Number 
formed from thickness of boundary layer and speed 
at the outer edge corresponds to the commonly used 
Reynolds Number for pipes formed from diameter and 
mean speed. We assume that there is a definite func- 
tional relation between the Reynolds Number for 
transition and the initial turbulence. 3 

SKIN FRICTION WITH EDDYING BOUNDARY LAYERS 

The skin friction on a plate when the flow in the 
boundary layer is eddying may be estimated by carry- 
ing over the results of measurements of skin friction 
in pipes (References 12 and 15). Experiments on 
eddying flow in pipes show that the force, F a , per 
unit area of wall is given by 



^ a = 0.045 P 



v V' 4 

UrJ 



(12) 



where p is the density, U the speed at the center of 
the pipe, v the kinematic viscosity, and r the radius 
of the pipe. It is assumed that the same equation 

a Schiller's experiments showed that the critical Reynolds Number was always the 
same for a turbulence produced in a given manner, although the measurement of 
turbulence was not quantitative. 



applies to skin friction on a plate with eddying flow- 
in the boundary layer, if U is interpreted as the speed 
at the outer edge of the layer and r is replaced by 5, 
the thickness of the boundary layer. 

In eddying flow in pipes the speed, u f varies across 
the cross section according to the law 



(13) 



where y is the distance from the wall of the pipe. 
This formula fits the experimental observations very 
closely except very near the center of the pipe and very 
near the wall. By assuming the same distribution of 
speed in an eddying boundary layer, replacing r by 5, 
and by using the Karman integral relation (5) together 
with the value of F a from (12) substituted for 

— n which applied only to laminar flow, the 

thickness, <5, of the boundary layer may be determined 
and F a evaluated in a more useful form. We find 



1 TT2 dS 1 TT AU 



and the integral relation becomes 



7 rp dS 23 T dU_ 



-0.045 



</ r V" 



The solution of this equation is 

5 5/4 [/ns^ = o.3536 r Ui j* i/ 37/T dx 

For skin friction on a plate £7= constant = U 0 

5 = O.37l(0 /6 at 1/5 

TP/ v V /6 
f.-0.0577.p-^ { m ) 



aoid 



(14) 



(15) 



(16) 



(17) 



It should be noted that F a , the force per unit area, is a 
function of x. We are more interested in the average 
force per unit area, F a> on a plate of length, I, namely 

i f F a dx which turns out to be 
I Jo 



f« = 0.072 P 



2 \Ul 



\i/ 6 

7 



The force coefficient, C F , is given by 

n F a 0.072 
CF ~i P U 2 ~ R 1/6 

where R is the Reynolds Number^ 



(18) 



(19) 
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The experiments of Wieselsberger (Cf. reference 15) 
on skin friction indicate good agreement with formula 
(19), except that a better value of the constant is 
0.074. In the calculations made in this paper on 
eddying boundary layers, we shall use the formula 



SEPARATION WITH EDDYING BOUNDARY LAYERS 

When the motion in the boundary layer becomes 
eddying, the phenomenon of separation is delayed. I 
In eddying motion, there is a more thorough mixing of 
the air particles and the driving action of the outer j 
layers is greater. This delayed separation produced 
by the eddying motion of the boundary layer is respon- 
sible for the great variation of the drag coefficient of 
spheres and cylinders in the critical region. 

DEFINITION OF TURBULENCE 

Before it is possible to speak of the measurement 
of turbulence in wind tunnels, we must be able to give 
a precise definition. The definition adopted is as 
follows: The turbulence at a given point is taken to 
be the ratio of the square root of the mean square of 
the deviations of the speed from its mean value to the 
mean value. At any point in the wind stream, the 
speed fluctuates with time about a mean value. The 
turbulence is the mean fluctuation taken in a definite 
manner and expressed as a percentage of the mean 
speed. A turbulence of 1 per cent means an equiva- 
lent sine wave fluctuation of ±1.4 per cent from the 
mean value. 4 

The adoption of this simple scalar quantity needs 
some justification since in theoretical treatments of 
eddying motion it is found necessary to separate the 
speed fluctuation into components. Experiment shows 
that the mean speed is very nearly constant over the 
greater part of the cross section of wind-tunnel air 
streams. The shape of the distribution curve changes 
as we go downstream only in the vicinity of the walls. 
There are accordingly no forces acting between adja- 
cent layers in the core of the air stream and we may 
assume that the fluctuations of speed are entirely 
random. In the neighborhood of the walls or close 
behind the honeycomb, this assumption is not true 
and the separate components of the velocity fluctua- 
tions as well as their phase relations must be considered. 
To characterize the air stream of a wind tunnel, we 
need only to consider the fluctuation of the absolute 
value of the speed so long as we do not get close enough 
to the honeycomb to be able to detect the honey com!) 
pattern in the distribution of mean speed. 

We need also to discuss the space distribution of the 
turbulence for many investigators have supposed that 

* This definition differs somewhat from that given in Reference 4. There, the 
double amplitude of the equivalent sine wave was given i. e. 2.8. Values in Refer- 
ence 4 should be divided by 2.8 to be comparable with those given in the present 
paper. 



some space characteristic of the turbulence was the 
important factor. It has been suggested for example 
that the turbulence should be characterized by the 
ratio of the diameter of an average eddy to the model 
diameter, a quantity related perhaps to the ratio of 
the diameter of a cell of the honeycomb to the model 
diameter. We find no evidence to support this idea. 
Experiment shows (reference 4) that the turbulence is 
very nearly constant across the cross section in the 
region where the mean speed is constant. It has further 
been shown (reference 16) that turbulence may be 
introduced without effect so long as the turbulent air 
does not reach the boundary layer at the surface of the 
model. Therefore unless the space distribution of the 
turbulence across the air stream is such that large 
changes take place in a distance comparable with the 
thickness of the boundary layer, we would not expect 
the space distribution to be an important factor. We 
expect the observed force to be governed by the turbu- 
lence of the air entering the boundary layer. 

The square root of the mean square deviation is 
chosen instead of the mean deviation for convenience, 
since the final measuring instrument is a hot wire 
alternating current milliammeter which gives this type 
of mean. 

DESCRIPTION OF WIND TUNNELS USED 

There are three wind tunnels at the Bureau of 
Standards, all of which were used in this investigation. 
Sketches of the three tunnels are shown in Figure 2. 
The following brief descriptions will serve to supple- 
ment the sketches. 

The 4%-foot tunnel is of the National Physical 
Laboratoiy type, octagonal in cross section, the 4K- 
foot dimension being between opposite faces. The 
faired entrance of the tunnel is about 4 feet long, the 
parallel portion 25 feet long and the exit cone 15 feet 
long. The diameter at the propeller end is 9 feet. 
The tunnel room is 68.5 feet long, 28.3 feet wide, and 
18 feet high, and is divided transversely near the exit 
end of the tunnel by a wall honeycomb, consisting of 
pasteboard tubes 1 inch in diameter and 4 inches long, 
packed in a light framework, which is covered on both 
sides by wire netting. This open honeycomb structure 
serves to damp out the swirl and eddies in the air 
stream as it returns from the fan to the tunnel entrance. 
A speed range of 25 to 110 feet per second is attained 
with the expenditure of from 2 to 75 horsepower. The 
propeller is 8 feet 11 inches in diameter and is driven 
at speeds from 170 to 870 revolutions per minute by a 
direct current shunt-wound motor. The tunnel is of 
wooden construction on a steel framework. 

The 3-foot tunnel is of the Venturi type, circular in 
cross section. The entrance cone is 12 feet long, the 
working portion 6 feet long, and the exit cone 33 feet 
long. The diameter at the end of the exit cone and at 
the front of the entrance cone is 7 feet. The room is 
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103.5 feet long, 21.3 feet wide and 14 feet high. The 
honeycomb at the tunnel entrance is of one-quarter- 
inch wood, the cells being 3 by 3 by 12 inches long. A 
plaster fairing is used in t lie entrance cone; otherwise the 
tunnel is of wooden construction throughout . The wall 
honeycomb is identical in type with that of the 4^-foot 
tunnel but is installed near the entrance end. Speeds 
up to about 240 feet per second are obtained with a 
motor rated at 100 horsepower. The propeller is 



(> 



current motor. The maximum speed obtained is 
approximately 100 feet per second. 

MEASUREMENT OF TURBULENCE 

The measurement of turbulence as defined in a pre- 
ceding section was made by the use of the hot wire 
anemometer in conjunction with an amplifier and 
apparatus for compensating for the lag of the hot wire. 
The early form and theory of the apparatus are given 
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Figure 2 



8-bladed, and is rotated at speeds up to 1,000 revolu- 
tions per minute. 

The 10-foot wind tunnel is outdoors and is of wooden 
construction. The cross section is circular and the 
length of the tunnel proper is 92 feet, A faired entrance 
bell, 4 feet long, is followed by a wooden honeycomb, 
cells 4 by 4 by 12 inches long. The cylindrical section 
is 50 feet long and the exit cone 34 feet long. The 
diameter at the exit end is 14 feet 2 inches. The 
4-bladed, 14-foot propeller is rotated at speeds up to 
550 revolutions per minute by a 200 horsepower direct- 
90351—30 2 



in reference 4. Since that report was prepared, sev- 
eral important modifications have been made in the 
interest of convenience, portability, and accuracy. 
These modifications are treated in an appendix to this 
paper, which in itself forms a supplement to reference 4 ■ 
A photograph of the amplifier and accessory apparatus 
in its modified form is given in Figure 3. 

The turbulence was measured at three sections of 
the 4K-foot tunnel as indicated in Figure 2 and at the 
working section of the 3-foot tunnel and 10-foot I unnrl 
The mean values obtained are given in Table I. 
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Table I 

Turbulence 

Location (per cent) 

3-foot tunnel, working section 0. 5 

10-foot tunnel, working section 1. 0 

4^-foot tunnel, downstream section 1. 2 

4}£-foot tunnel, working section 1. 6 

4H-foot tunnel, upstream section 2. 3 

Many repeat measurements show that the precision 
of the above values is of the order of 0.1 to 0.2. It is 



We have made tests on three spheres, of diameters 4, 
5, and 8.6 inches, at the five locations at which the 
turbulence was measured, except that measurements 
were not made on the large sphere in the 3-foot wind 
tunnel because of the large ratio of the diameter of the 
sphere to the diameter of the tunnel. Sketches of the 
suspensions used are given in Figure 4 and the results 
of the measurements are given in Figures 5 and 6. 




Figure 3 



obvious that the effects of turbulence could be studied 
over a reasonable range of values. 

MEASUREMENTS ON SPHERES 

The effect of turbulence in wind-tunnel measure- 
ments was first discovered for spheres and it has often 
been suggested that measurements on a sphere be 
used as a measure of turbulence. O. Flachsbart (refer- 
ence 17) in a discussion of this proposal points out the 
necessity of some standard method of suspension, if 
comparable results are to be obtained. We have not 
used exactly the same method of measurement through- 
out our sphere tests but Flachsbart's results for the 
suspensions used in our tests show that the differences 
introduced by this fact are small. 



The 4-inch sphere was mounted on two wires ar- 
ranged in the form of a V in a plane transverse to the 
wind direction. The force was computed from the 
deflection of the sphere as a pendulum and the small 
wire corrections were determined by computation. 
Some experiments were made on the 5-inch sphere in 
the 4%-foot tunnel on a bell crank. Correction was 
made for the motion of the scale pan and the spindle 
correction was determined with the sphere mounted 
in front of the spindle, We prefer as a standard 
method the arrangement used for the 5-inch sphere in 
the 3-foot and 10-foot tunnels and for the 8.6 sphere in 
all cases, namely, a tail spindle supported by two 
V's with a shielded counterweight. In this arrange- 
ment all wires are behind the sphere. The force is 
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computed from the downstream deflection of the sus- 
pended system. The spindle correction is determined 
in the same manner, the sphere being detached and 
supported in front of the spindle. A comparison of 
Figures 5 and 6 shows a difference in shape which we 
have traced to an effect of the balance windshield. 

It will be noted that the curves run from right to left 
in order of increasing turbulence in accordance with 
the interpretation of Wieselsberger's measurements 
(reference 18) of the effect of screens. To use the 




Mounting used on 6.6" 
sphere throughout and 
on 5" sphere in wind 
tunnels E and 3 




Mounting used 
on 4" sphere 




Pivot 



Top of tunnel 
Mounting used on 5" sphere in Wind tunnel I 

Wind shield for arm not shown > 



O 



Figure 4— Mountings used for the measurement of the air resistance of spheres 

sphere as an instrument for the quantitative measure- 
ment of turbulence, we must give some more precise 
definition of the critical Reynolds Number than has 
hitherto been given. We suggest that the critical 
Reynolds Number be defined for this purpose as the 
value of the Reynolds Number at which the drag 
coefficient of the sphere is 0.3. For a given condition 
of turbulence the results for the different spheres give 
slightly different values of the critical Reynolds Num- 
ber traceable in part to the differences in mounting, 
but the extreme difference is only 8 per cent and the 
mean difference much less. Adopting the mean values 
we obtain the values given in Table II. 



Table II 

Turbulence Critical Reynolds 

(per cent) Number for sphere 

0.5 270, 000 

1.0 232,000 

1.2 197,000 

1.6 164,000 

2.3 132,000 

These results are plotted in Figure 7, which is a cali- 
bration of the sphere as an instrument for measuring 
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Figure 5.— Drag coefficients for 4" and 5" spheres. p = air density, t/ 0 =air 
speed, ^L=area of cross section of sphere = 7r~— D=diameter of sphere, v = 
kinematic viscosity 

turbulence. The circles show the approximate limits 
of individual values. The accuracy while not high is 
perhaps sufficient for most purposes. 
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Figure 6— Drag coefficients for 8.6" sphere, p=air density, I7.=air speed, 
yl=area of cross section of sphe^e= 1 ^- , Z>=diameter of sphere, f = kinematic 
viscosity 

DISCUSSION OF SPHERE MEASUREMENTS 

The interpretation of the sphere measurements in 
the light of the boundary layer theory is as follows. 
At low values of the Reynolds Number of the sphere, 
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the flow in the boundary layer is laminar and separa- 
tion takes place in a manner governed by the laws of 
laminar flow. As the Reynolds Number of the sphere 
is increased, the Reynolds Number formed from the 
thickness of the boundary layer and the speed at the 
outer edge reaches the value at which eddying How be- 
gins at some point upstream from the separation point. 
Separation is then delayed, the wake is smaller and the 
form resistance is decreased. The skin friction is such 
a small part of the total drag that the effect of the in- 
creased skin friction due to the change from laminar to 
eddying flow is inappreciable. 

If the air stream is initially turbulent the change will 
occur at a lower value of the Reynolds Number of the 
boundary layer and hence the critical Reynolds Num- 
ber for the sphere will be reduced by an amount which 
increases with increasing turbulence. 

We can not make the interpretation a quantitative 
one until more satisfactory methods are developed for 
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Figure 7.— Critical Reynolds Number of spheres (at which C A (figs. 5 and 
6) is 0.3) as a function of the turbulence 

treating the phenomenon of separation with eddying 
as well as laminar boundary layers. An approximate 
semiempif i< al treatment is given by Ono (reference 19) 
but Tollmien (reference 20) questions the legitimacy 
of t he approximations used. 

MEASUREMENTS ON STREAMLINE BODIES 

Having obtained a very good correlation between 
measurements on spheres and the turbulence as meas- 
ured by our hot wire anemometer and associated appa- 
ratus, experiments were begun on streamline bodies. 
The first model used was a bomb model known as 
1I-Q-12 which we had at hand. We were astonished 
to find no large effect in view of the results on the 
N. P. L. models in different wind tunnels, for we had 
expected that generally similar results would be ob- 
tained on all streamline bodies. The actual bomb 
model departed in several apparently minor respects 
from a good streamline body and we therefore made a 
wooden replica without protuberances of any kind. 
The shape of this replica is shown in Figure 8 and the 
results of measurements of the air force at the five 
measuring stations are given in Figure 9. These 



curves show certain systematic changes but practically 
all of the points are within 10 per cent of a mean curve. 

The forces were measured by swinging the model on 
four wires arranged in two V's from the nose and tail, 
measuring the deflection at several wind speeds, the 
weight, and wire lengths, and computing the total drag 
of model and wires. The drag of the wires, amounting 
to about 75 per cent of the drag of the model, was 
computed, since we have found this procedure to give 
accurate results. To minimize errors due to the cor- 
rection for the wires, we have used the same wires at 
all stations in the 4%-foot tunnel and in the 10-foot 
tunnel. In the 10-foot tunnel, the model was hung 
from an auxiliary frame supported in the wind stream. 

0.0125 ft. ■ 



.469 ft.^ 




1 4.972 ft. 




H-Q-12 




Area of cross -section 


0.0873 ft* 


Length 




1.6470 ft. 


Mox. diameter 


0.3330 ft. 




Figure 8. 



N.P.L. shorf model 
Votum e 0. 0385 ft. 3 

(Volume) 2 ' 3 0.2677 ft. 2 

Length 2.3340 ft. 

Mox. diameter O.3500 ft. 



N.P.L. long model 
Volume 0. 1886 ft. 3 

(Volume)*' 3 0.3291 ft. 2 

Length 2.8600 ft. 

Mox. diameter 0.3500 ft. 

-Longitudinal sections of bodies of revolution on which force measure- 
ments were made 



Care must be taken in the interpretation of curves 
such as are shown in Figure 9. The observations are 
made under such conditions that the forces and speeds 
are subject to errors of roughly constant absolute 
magnitude. The percentage error is therefore much 
greater at the lower speeds. Under these conditions 
a plot of coefficients may prove misleading unless one 
remembers continually that the experimental errors 
are different in different parts of the diagram. We 
estimate that the probable errors range from some 
13 per cent at a speed of 20 ft./sec. to perhaps 2 or 3 
per cent at a speed of 100 ft./sec. 

No corrections have been applied for pressure drop, 
because we do not believe that the method of cor- 
rection is as yet well established. We give neverthe- 
less the data required for making the correction. 
Table III gives the mean decrease in static pressure per 
foot length divided by the velocity pressure for the 
several stations. 
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Table III 

Station Pressure drop 

3-foot 0. 007 

10-foot • 002 

4^-foot, downstream section . 0035 

4}Hoot, working section , 0075 

4}£-foot, upstream section . 0035 

The volume of II-Q-12 is approximately 0.077 
cubic foot and the area of cross section, which is used 
in defining the drag coefficient for this model, is 0.0872 
square foot. Hence the corrections for pressure drop 
applied in the usual manner to the coefficients of 
Figure 9 are to be found by multiplying the values 

0 077 

of the pressure drop given in Table III by or 

0.885, giving 0.006, 0.002, 0.003, 0.007, 0.003, or 
approximately 10, 3, 5, 12, and 5 per cent for the 
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Figure 9.— Drag coefficients for II-Q-12 wooden replica. P =air density, = Uo 
air speed, vl=area of cross section (maximum section), Reynolds Number= 

— where / is the length of model, and v is the kinematic viscosity 

V 

several stations. The application of these corrections 
would not change the conclusion that the effect of 
turbulence is small compared with the magnitude of 
the discrepancies found between the results on the 
N. P. L. models in the several wind tunnels. 

It was obviously desirable to pass immediately to 
models which gave a larger drag coefficient in some 
wind tunnels than in others and we accordingly con- 
structed a wooden replica of the short airship model 
used in the N. P. L. comparative tests. The results 
of a large number of runs are shown in Figure 10. 
Although it is difficult to accurately measure the small 
forces, there is clearly a very large difference between 
the measurements at the upstream and downstream 
90351—30 3 



sections of the 54-inch wind tunnel. The variation 
may be attributable in part to a change in the charac- 
ter of the surface of the model or to a change in shape 
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Figure 10.— Drag coefficients for wooden replica of N. P. L. short model, p=air 
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of the model inasmuch as the tests covered a consider- 
able period of time. To exclude this complication so 
far as possible, arrangements were made to use the 
U. S. Navy replicas of the N. P. L. models. The 
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Figure 11 —Drag coefficients for metal replica of N. P. L. short model, p«air 

density, U 0 =air speed, Reynolds Number=^ where I is the length of the 

model, and v is the kinematic viscosity 

nominal dimensions of the models are given in Table IV 
and the form is shown in Figure 8. The results of 
measurements of the drag are shown in Figures 11 
and 12. 
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TABLE IV.— NOMINAL DIMENSIONS OF U. S. N. 
REPLICAS OF N. P. L. MODELS 



Distance 
along axis 
from nose 


Short model 
moan radius 


Long model 
mean radius 


0 


0 


0 


0.250 
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We may sum up the experimental results on the. 
effect of turbulence on the resistance" of streamline 
models by saying that some models show small effects, 
others show large effects, and that the large effects are 
confined to a certain range of values of the turbulence. 
For example, the N. P. L. models show small effects 
if the turbulence is less than about 1.3 per cent, large 
effects if the turbulence is greater, the effect of in- 
creasing turbulence being to increase the insistence. 

We have then an explanation of the results of the 
American tests on these models. The old variable 



density tunnel was very turbulent as indicated by the 
low value of the critical Reynolds Number for spheres 
(about 94,000, Cf. reference 21). In agreement with 
this indication, the measured drag of the airship models 
was greatest in the variable density tunnel. (Refer- 
ence 3.) At the other extreme, the lowest values of 
the measured drag of the airship models were obtained 
in the 3-foot and 10-foot wind tunnels of the Bureau of 
Standards, which have small turbulence. In the new 
variable density wind tunnel the turbulence has been 
greatly decreased (reference 23) and it would be ex- 
ceedingly interesting to determine whether the airship 
models now give much lower drag coefficients. 

DISCUSSION OF MEASUREMENTS ON STREAMLINE 
BODIES 

The experimental results on airship models leave one 
in a very confused state of mind as to the interpreta- 
tion of model experiments on airship hulls and as to 
the explanation of the puzzling feature that the effects 
are large only for some models over a certain range of 
ralues of the turbulence. We do not desire to leave 
the subject in this state and we believe that the 
boundary layer theory as we have outlined it will 
account for the observed facts. Because we are 
handicapped by the absence of methods of exact 
mathematical treatment, we can only give a discussion 
based on rather crude mathematical approximations 
with the hope that the reader will obtain some sug- 
gestion as to how the observed results may follow from 
the boundary layer theory. 

We first determined experimentally that the phenom- 
enon of separation does not enter. This was done by 
placing a thin film of oil on the surface and noting that 
there was no region of reversed flow. We therefore 
expect that the form resistence will be small. In fact, 
pressures have been measured on a model which is 
substantially the N. P. L. long model (reference 22) 
and it was found that the form resistance was prac- 
tically zero. In other words, the observed resistance 
is due almost entirely to skin friction. We shall adopt 
this assumption and attempt to calculate only the skin 
friction. 

We propose to compute the skin friction of the 
"equivalent flat plate/' i. e. on a section of a two- 
dimensional plane surface, the width of the section at 
a given distance from the front edge being equal to 
the circumference of the model at the same distance 
from the nose of the model, and the speed at a given 
distance from the front edge of the plate being the 
same as the speed computed by Bernoulli's theorem 
from the pressure distribution at the same distance 
from fche nose of the model. We shall identify the 
skin friction on the equivalent plate with the resistance 
of the model. The following assumptions are implied 
in the above procedure. 

1. No distinction is made between distances along 
the surface and distances along the axis. 
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2. The cosine of the angle between any surface ele- 
ment and the axis is taken equal to unity. 

3. The thickness of the boundary layer at a given 
point is assumed small compared with the radius of 
cross section of the model at the same point (obviously 
not true far back on the tail). 

4. The equations of two-dimensional flow are used. 
For the computation of the point at which the flow 

in the boundary Layer changes from laminar to eddy- 
ing, and for the computation of the skin friction on 
the part of the surface for which the flow in the 
boundary layer is laminar, w r e use the formulas of 
equations 7 and 9. These are based on the use of a 
linear velocity distribution in the boundary layer and 
the integral relation of Karman. 

For the computation of the skin friction on the part 
of the surface for which the flow in the boundary layer 
is eddying, we use equation 20 together with an assump- 
tion made by Prandtl (Reference 15) in connection with 
skin friction on plates, namely, that the force is the 
same as if the flow were eddying from the beginning. 
If the flow were eddying from the beginning, the co- 
efficient c 2 , for a surface of length l 2 in a stream of 
speed U 0 , would be given by 



c 2 = 



where R 2 = ^^* 



For a surface of length l u the coefficient, c x is given by 
0.074 , „ Uoh 



Ci = 



R^ 



where R\ •■ 



Hence in the first case, we may find the desired average 
coefficient, c, for that part of the surface between l x 
and l 2 (both k and l 2 being measured from the nose) 
by stating that the weighted average of c x and c must 
be c 2 , or 

Cili±c(l 2 — 1\) =*c 2 l 2 



whence 



C 2 l 2 — C\l\ 



(21) 



We add another to our formidable list of assumptions 
and approximations by not considering the variations 
of the speed at the outer edge of the boundary layer 
in this computation, taking the wind tunnel speed, U 0 , 
to compute c 2 . 

We are now prepared to compute the skin friction 
on the N. P. L. long model for wind streams of differ- 
ing amounts of turbulence in accordance with the ideas 
outlined in the section on eddying boundary layers, 
namely, that to each degree of turbulence there corre- 
sponds a different value of the critical Reynolds Num- 
ber of the boundary layer for which the flow changes 

JJ 

from laminar to eddying. The values of jj computed 



by Bernoulli's theorem from the pressure distribution 
given in Reference 22 are shown in Figure 13. We 
prefer to write equation (7) as 



where 



(22) 



*=7A4X0 d *- (23) 

Wo/ 

We note that the Reynolds Number of the boundary 
layer, R s , is 



1.0 



.4 
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Figure 13.—^ and(~^y for N. P. L. long model computed from pressure 
distribut ion of Reference 22 

Introducing R= J ~^j the Reynolds Number 8 of the 
model, where / is the length of the model, wc find 
Rs U /12~ 



\R 



K. 



(25) 



jj- ) vs x. The value of K 

was found by graphical integration as a function of x } 

R 

and from the values of K , as a function of x } which 

is plotted in Figure 14. Equation (9) for the laminar 
force coefficient may be written 

U 



. xT dX 



r (1 



(26) 



We shall eventually need the contribution to the force 
coefficient, C v of the model defined in terms of (vol- 
ume) % , for which purpose C F must be multiplied by 
the surface area, A Xy of the model from the nose to the 

« The reader should note the use of the length of the model instead of (volume)"^ 
in the definition of R. It will be appreciated that the length gives a better basis of 
comparison. 
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point x, and divided by (volume) * A . As before, we 
introduce the Keynolds Number of the general flow 
and write 

Cr^R A z _ II A x C*Uo^ (27) 
(Vol) « V 3 x J„ JK * 
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Figure L4.— and ( -^yJ- 3 V« for N. P. L. long model. « a is Reynolds 

Number of boundary layer at a distance x from the nose, R is the Rey- 
nolds Number of the body, CV is the average laminar coefficient of fric- 
tion, ylx is the surface area of the model from the nose to the point x 



This quantity, which is to be divided by ^R to give 
the laminar contribution to the total force coefficient, 
is plotted as a function of x in Figure 14. 

We next compute a table showing the laminar con- 
tribution to the total force coefficient for varying R 
and x, a table showing the contribution of the surface 
between x and I for which the flow is eddying (com- 
puted from Equation 21 by multiplying c by the sur- 
face area and dividing by (volume) for varying 
R and x } and finally a table of computed total force 
coefficients (adding corresponding entries in the first 
two tables) for varying R and x. Parts of these tables 
for the N. P. L. long model are given in Tables V, VI, 
and VII. Table VII permits the easy calculation of 
the total drag coefficient, when the point at which the 
flow changes from laminar to eddying is known. If, 
for example, the change always occurs at a fixed point, 
values are taken from the horizontal line corresponding 
to that point. 

Assuming the change from laminar to eddying flow 
to occur when the Reynolds Number, R Sy of the bound- 
ary layer is equal to 1,500, 2,000, 2,500, 3,000, and 
3,500 we find x from Figure 14 and compute the curves 
shown in Figure 15. The computations for R 5 equal 
to 1,500 and 3,500 are shown in part in Table VIII. 




5 6 7 8 9/0 6 



7 8 9I0 6 2 3 4 5 6 7 8 9I0 7 

Reyno/ds Number, H 

Figure 15.— Computed drag coefficients for N. P. L. long model for flow changing from laminar to eddying under different conditions. Upper line is for eddying 
flow throughout, lower line for laminar flow throughout. The dotted lines are for transition at the indicated distances from the nose. The curves are for transi- 



tion at the values of the Reynolds Number, Bg of the boundary layer indicated C,— 
the length of the model and v the kinematic viscosity 



Force 
&p(VoI)*/» C/.i 



p=air density, l/*=air speed. R=*~ * where / is 
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TABLE V.— LAMINAR CONTRIBUTION TO TOTAL FORCE COEFFICIENT, N. P. L. LONG MODEL 



Speed, Uo, ft./see 

Reynolds Number, RxlO- 3 . 



1. 10, 
.15, 
.20, 
.30, 
.40, 
.50, 



10 
180 



I. OIL'L' 
.0031 
.0039 

. 0055 
.0070 
. oosi 
.0097 
.0117 
.0134 
. 0150 
.0101 
.0177 
.0188 
. 0190 

. 0220 



0. 0015 
.0022 
.0028 
. 0030 
.0050 
.0060 
. ()069 
.0083 
.0095 

.0100 
.0115 
. 0125 
.0133 
.0141 
. 015(5 



540 



0. 0012 
.0018 
.0022 
.0032 
.0041 
.0049 
. oor.ii 

.OOOS 

. 007 s 
.0086 
. (10!).', 

.0102 
.0109 
.0115 
.0127 



40 
720 



0.0011 
.0015 
. 0019 
.0(12* 
. 0035 
.0012 
.0048 
. 0059 
.0007 
. 0075 
. 00S2 
. (Miss 
.0094 
.0100 
.0110 



50 

•.ton 



0.0010 
.0014 
. 0017 
.0025 
. 0031 
.0038 
.0013 
.0053 
. 0060 
.0067 
.0073 
. 007! I 
.0084 
. OUST 
. 00! »s 



60 
1080 



0. 000!) 

.0013 
. 00 If. 
.0023 
. 0029 
.0035 

. 0010 

.0048 
.0055 
.0061 
.0067 
.0072 
.0077 
. OOM 
. 0090 



70 
1260 



o. onus 
.0012 
. 0015 
.0021 
. 0027 
.0032 
. 0037 
.0044 
.0051 
.0057 
.0002 
.0067 
.0071 
.0075 
. imk; 



80 

HID 



0. OOOS 

.0011 
.0014 
.0020 
. 0025 
.0030 
.0034 
.0042 
.00 is 
.0053 
. oor.s 
.0063 
. 0067 
.0071 
.007s 



90 
1620 



100 
1800 



0. 0007 
.0010 
. 0013 

. 0010 

.0023 
.0028 
. 0032 
.003!) 

.ooi:. 
.0050 
.0055 
.0059 

.0002 
.0067 
.0073 



0. 0007 

.0010 

.0012 
.0018 
. 0022 
.0027 
.0031 
.0037 
.0043 
.0047 
. 0052 
.0056 
.0060 
.0063 
.0070 



TABLE VI. — EDDYING CONTRIBUTION TO TOTAL FORCE COEFFICIENT, N. P. L. LONG MODEL 



x ft. 



Speed, Uo, ft./sec 

Reynolds Number, RxlO" 3 . 



1. 10, 
.15. 
.2__ 
.3.- 
.4_. 
.5__ 
.6.. 



10 
180 



0.0472 

.0100 
.0450 
.012s 
. 0406 
.0383 
. 03;,!) 
.0316 
.0273 
.0231 
.0191 
. 0152 
.0113 
.007!) 
.0010 



20 
360 



0.0411 

.0100 

.0391 
.0372 
.0353 
.0333 
.0312 
.0274 
.0238 
. 0201 
. 0166 
.0132 
. OOOS 
. 000!) 
. 1100!) 



30 
540 



0. 0378 
.0369 
. 0360 
. 0342 
.0325 
. 030(5 
.02 s 7 
.0253 
.0219 
.0185 
. 0153 
.0121 

. 00!) 1 

.0064 

.OOOS 



40 

720 



0.03.'.!) 
. 0350 
. 0312 
. 0326 
. 0309 
.0291 
. 0273 
.0240 
. 0207 
.0176 
.0145 
.0110 
. OOSO 
. 0000 

.0007 



50 
900 



0. 0342 
. 0331 
. 032(5 
.0311 
.0294 
.0278 
.0260 
.022!) 
.0198 
.0167 
.0138 
.0110 

. 00S.3 
. 00 5S 
.0007 



60 
1080 



0. 0321) 
.0320 
.0312 
.0298 
. 02S2 
. 02(5(5 
.0250 
.0219 

.011)0 

.0160 
.0133 
.0105 
.0079 
.0055 
.0007 



70 
1260 



0. 0319 
.0311 
. 0301 
.0289 
.0274 
.0259 
.0243 
.0213 
.0185 
.0150 
.0129 
.0102 
.0077 
.0054 
.0006 



1440 



0. 0309 

. 0302 
. 0204 
.0281 
. 0265 
.0251 
.0236 
.0207 
.0179 
.0152 
.0125 
.0100 
.0074 
.0052 

.0000 



90 
1620 



0.0304 
.0297 
.0290 
.0276 
.0262 
. 0217 
.0233 
. 0203 
.0176 
. 01 10 

.0123 
.0098 
.0073 
.0051 
. oooo 



100 
1800 



0. 0298 
.0290 
. 02S3 

.0270 
.0256 
. 0212 
.0227 
.0199 
.0172 
. 0146 
.0120 
.0095 
.0072 
.0050 
.0006 



TABLE VL-TOTAL FORCE COEFFICIENT FOR^L W^MODE L, W,TH FLOW CHANGING FROM 



X ft. 



Speed, Uo, ft./sec 

Reynolds Number, Rxl0-3_ 



1 10. 
.15. 
.2.. 



10 
180 



0.049! 
.0491 
.0489 
.0483 
.0476 
. 0107 
.0456 
.0433 
.0407 
.0380 
.0355 
.032!) 
.0301 
.0278 
.0230 



20 
360 



0. 0120 
.0422 
.0419 
.0411 
. 0103 
.0303 
. 03S1 
.0351 
.0333 
.0307 
.0282 
.0257 
.0231 
. 0210 
.0165 



30 
540 



0. 0390 

.0387 
.0382 
.0371 
. 0366 
.0355 
.0343 
. 0321 
.0297 
.0271 
. 021 s 
.0223 
.0199 
.0179 
.0135 



40 
720 



0. 0370 
. 0305 
. 0361 
.0354 
.0344 
.0333 
.0321 
.0299 
.0274 
.0251 
.0227 
.0204 
.0180 
. 0159 
.0117 



50 
900 



0. 0352 
.03 IS 
. 0343 
.0336 
.0325 
.0317 
.0303 
.0282 
.0258 
.0234 
.0211 
.0189 
.0167 
.0147 
.0105 



60 
10M> 



0. 033 s 
. 0333 
. 03 2N 
.0321 
.0311 
.0301 
.021)0 
.0207 
.0245 
.0221 
.0200 
.0177 
.0150 
.0136 
.0097 



70 


80 


90 


100 


1260 


1440 


1620 


1800 


0. 0327 


0. 0317 


0.0311 


0. 0305 


.0323 


.0313 


.0307 


. 0300 


.0310 


.0308 


.0303 


.0295 


.0310 


.0301 


. 0205 


.0288 


.0301 


. 02D0 


. 02S5 


. 027s 


.0201 


.0281 


.0275 


. 0200 


. 02S0 


.0270 


.0265 


. 02 66 


.0257 


. 0240 


.0242 


.0236 


.0236 


.0227 


.0221 


.0215 


.0213 


.0205 


.0100 


.0193 


.0190 


.0183 


.0178 


.0172 


. 0169 


.0163 


. 0157 


.0151 


.0148 


.0141 


.0135 


.0132 


.0120 
.0089 


.0123 


.0118 


.0113 


.0084 


.0079 


| .0076 



TABLE Vill.-COMPUTED DRAG COEFFICIENTS, 
N P L LONG MODEL, WITH CHANGE FROM LAMI- 
NAR TO EDDYING 'FLOW IN THE BOUNDARY 
LAYER DETERMINED BY Rs=l,o00 AND 3,500 



R 


tf 5 = 1,500 


i? 5 =3,500 


** 

4r 


(Fig* 14) 


C, 

(Table 7) 


■ylR 


(Fig 2 14) 


c, 

(Table 7) 


90,000 

ISO. 000 
360,000 
540.000 
720,000 

000.000 

1,080,000 
1,260,000 
1,440,000 
1,620,000 
1, 800, 000 


5.00 
3. 53 
2.50 
2.04 
1.76 
1.58 
1.44 
1.34 
1.25 
1. 18 
1. 12 


2.60 
2.00 
.95 
.75 
.61 
.52 
.44 
.40 
.36 
.33 
.31 


0. 0321 
.0278 
. 0337 
.0326 
.0320 
.0314 
. 0307 
.0301 
. 0204 
.0292 
.0287 


11.65 
8.25 
5. 82 
4. 76 
4. 12 
3.68 
3.37 
3.12 
2.91 
2. 75 
2.61 


2.53 
2. 2s 
2.08 
1.89 
1.50 
1.20 
1.07 
1.00 


0.0313 
.0221 
.0157 
. 0139 
. 0137 
.0141 
.0147 

.oiNo 

.0205 
.0214 
.0215 



Figure 15 gives a general survey of the computed 
coefficients over a wide range of Reynolds Number. 
The curves (straight lines on the logarithmic plot) for 
completely laminar and completely eddying flow are 
shown as the limiting cases. In addition two dotted 
lines are drawn to show the curves that one obtains 
for the flow changing at a definite position on the 
model. The part of the curves for Reynolds Numbers 
up to 2.5 x 10 6 that we have covered in our experi- 
ments, is replotted in Figure 16 on a nonlogarithmic 
plot for direct comparison with the experimental 
results given in Figure 12. 

A comparison of Figures 12 and 16 shows a general 
agreement in order of magnitude of the forces and in 
shape of the curves, except for the three lower curves 
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of Figure 12 below 50 ft. /sec. In order to account for 
the observed shape of these curves and in view of other 
evidence that will be presented later, we have been led 
to the assumption that the flow soon becomes eddying 
behind the maximum cross section or speaking more 
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Figure 16.— Reproduction of a part of Figure 15 on a nonlogarithmic scale with 
additional curves for transition from laminar to eddying flow at specified dis- 
tances from the nose 

exactly that the contraction of cross section introduces 
turbulence, producing the same effect as any upstream 
turbulence on the flow. We do not think it unreason- 
able to suppose that the necessary contraction in cir- 
cumference and increased thickening of the 
boundary layer resulting from the reduction 
in cross section of the body does not take 
place smoothly and uniformly; that there 
is a folding or wrinkling of the layer which 
produces disturbances of the same nature 
as the turbulence of the wind tunnel air 
stream. Let us trace through on Figure 
15 the consequences of such an assumption. 
We suppose for simplicity that at the point 
x =1.9 feet which is some distance behind 
the maximum cross section, turbulence 
arises which would cause a change from 
laminar to eddying flow at a Reynolds Num- 
ber of the boundary layer, R h of 2,000. 
We suppose further that the turbulence of 
the wind tunnel air stream corresponds to 
an R s of 3,500. Beginning at low values 
of the Reynolds Number of the body, i. e. 
at the extreme left, we have laminar flow 
throughout and follow the line of laminar flow. As 
soon as we strike the curve for R 6 = 2,000, the flow 
becomes eddying on a part of the tail and we follow 
the curve for R b = 2,000 (the point of transition moving 
forward) until we reach the line corresponding to the 
point of transition, 35=1.9, the lower of the dotted lines. 
The point of transition then remains stationary, and 
we follow the dotted line for x = 1.9 until R s at this 



point reaches 3,500 as indicated by intersection with 
the curve for R 6 = 3,500. We then follow the curve 
for R$ = 3,500. A consideration of Figure 16 on which 
several curves for different points of transition are 
dotted in will show how the experimental curves can 
be fitted by a proper choice of point of transition curve 
and #5 = constant curve. The transition from one 
curve to another is undoubtedly not as sharp as indi- 
cated by this simple analysis. 

There is one other feature of the curves of Figure 12 
not yet accounted for, namely, the fact that for turbu- 
lence below about 1.3 per cent, there is apparently no 
effect of turbulence. We believed at first that this was 
caused by the disturbance from the supporting wires 
at the nose introducing a turbulence of approximately 
1.3 per cent, so that although the turbulence of the 
air stream was reduced below this amount, the bound- 
ary layer was always subject to the approximately 
constant turbulence introduced by the wires. If this 
interpretation is correct we should be able to obtain 
lower coefficients in the 3-foot wind tunnel by elimi- 
nating the disturbance at the nose. Figure 17 shows 
the results obtained by mounting the wooden replica 
of the N. P. L. short model on a tail spindle similar 
to the sphere mounting (fig. 4), compared with results 
on the same model and tail spindle with the front V 
at the nose. The effect is not large although the results 
with the V at the nose are somewhat higher and scatter 
more than those with all wires behind the model. 
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Figure 17.— Effect of nose suspension of the wooden replica of the N. P. L. short model and tail 
spindle on the drag coefficient. The values are higher than those of Figure 10 because no correc 
tion is made for the effect of the tail spindle 

The other possible explanations are of the same 
nature, for example, that turbulence is set up 
by the fore-and-aft oscillation of the model, or as 
van der Hegge-Zijnen has suggested (reference 11) 
that some turbulence is set up by the form of the 
nose. We have not as yet examined these possi- 
bilities, and must leave the matter open at the present 
time. 
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Calculations similar to those outlined for the N. P. L. 
long model were made for the short model and for 
II-Q-12, the pressure distribution being obtained 
experimentally in the 3-foot wind tunnel. It is 
unnecessary to give these computations in detail. It 
was found that curves very similar to those of Figure 
15 with just as much spread were obtained. This 
result, not in accordance with the experimental results 
on II-Q-12 (fig. 9), showed that the nature of the 
pressure distribution could not account for the small 
effect of turbulence on II-Q-12. We were forced to 
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Figure 18.— Drag coefficient for wooden replica of N. P. L. short model 
with wire rings 0.23 foot aft of nose. Figures on the curves give the 
diameter of the wire used to make the ring 

some other explanation and were led finally to the 
hypothesis already referred to, that the turbulence was 
introduced by the diminishing cross section of the body. 
The distinguishing feature of II-Q-12 causing it to be 
much less sensitive to changes of turbulence was seen 
to be the very forward position of the maximum cross 
section. Referring to Figure 16, for example, if at 
z = 0.5, corresponding to the upper dotted line, turbu- 
lence is introduced giving an jR 5 of 1,500, all of the 
observations must lie between the curve # 5 = 1,500, the 
dotted line, and the upper line for completely eddying 
flow. 

To obtain some further evidence that the diminution 
of the cross section introduced turbulence, experiments 
were made in which wire rings were placed around the 
model at several positions. The theory of these experi- 
ments is that if the wire is placed in a region where the 
flow is already eddying, the effect on the resistance of 
the model will be small, whereas if the wire is placed 
in a region where the flow is laminar, turbulence will 



be introduced and the resistance will be sensibly in- 
creased. Wires of different diameters produce differ- 
ent degrees of turbulence and even when the wire is 
placed in a region of eddying flow, there will be a slow 
increase of resistance with increasing diameter of the 
wire due to the resistance of the wire itself. Returning 
to Figure 15, suppose a wire placed at a; = 0.5 foot of 
such diameter as to produce a turbulence corresponding 
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Figure 19.— Drag coefficient for wooden replica of N. P. L. short model 
with wire rings 1.12 feet aft of nose 

to # 6 = 2,500. Suppose further that the turbulence of 
the wind tunnel corresponds to # 5 = 3,500 and that at 
x ==1.9 feet, owing to the diminishing cross section, a 
turbulence corresponding to #$=1,500 is introduced. 
The resistance coefficient would follow the curve R* = 
1,500 to its intersection with the curve for transition at 
1.9, follow the latter to its intersection with #5 = 2,500, 
follow #6 = 2,500 to its intersection with the curve for 
transition at a; = 0.5, follow this curve to its intersec- 
tion with # 5 = 3,500, and finally follow E 5 = 3,500. The 
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Figure 20.— Drag coefficient for wooden replica of N. P. L short model 
with wire rings 1.63 feet aft of nose 

Reynolds Numbers covered by our experiments do not 
permit the tracing of the last two stages. 

Some results for the wooden replica of the N. P. L. 
short model are given in Figures 18, 19, 20, and 21. 
Figure 18 shows that rings on the nose give curves 
similar to those obtained in wind tunnels of different 
degrees of turbulence. Figures 19 and 20 show that 
wire rings behind the maximum cross section give rela- 
tively small effects. The cross plot of coefficient 
against wire diameter for a speed of 80 ft. /sec. shown 
in Figure 21 shows the matter somewhat more clearly. 
The sharp rise in the curve for J- = 0.23 foot is inter- 
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preted as due to the effect of the turbulence set up by 
the wire. The slower rise is attributed to the increas- 
ing resistance of the wire itself. 

Figures 22, 23, and 24 give the results of similar 
measurements on II-Q-12. Here also, the effect of 
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Figure 21.— Drag coefficient for wooden replica of N. P. L. short model 
at an air speed of 80 ft./sec. with wire rings of varying diameter and 
position on the model. x= distance of ring aft of nose 

wire rings behind the maximum diameter is small and 
the effect of nose rings is similar to the effect of in- 
creasing wind tunnel turbulence. 

The use of wire rings has been suggested for routine 
measurements for the purpose of producing eddying 
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Figure 22.— Drag coefficient for wooden replica of II-Q-12 model with 
wire rings 0.164 foot aft of nose 

flow such as exists at high Reynolds Numbers. We 
may state the following conclusions as to this pro- 
cedure. A wire approximately 0.015 inch in diameter 
is required to give the full effect of a very turbulent 
air stream. The wire should be placed well forward 
on the model, but probably not so far forward as to be 



in the region of reduced speed. (Not less than 0.2 foot 
from the nose on the N. P. L. long model, for example. 
See Figure 13.) The wire should have relatively little 
effect in an air stream that is already turbulent. 

In closing this discussion of results on streamline 
bodies, we wish to exhibit a set of curves, Figures 25, 
26, and 27, computed by the methods previously out- 
lined for the three stations in the 4%-foot wind tunnel 

./Or 
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Figure 23.— Drag coefficient for wooden replica of II-Q-12 model with 
wire rings 0.97 foot aft of nose 

in comparison with the experimentally measured curves. 
We have adopted for each model a point at which we 
assume turbulence to be introduced by the diminish- 
ing cross section, and for each station a value of R 6 
characteristic of the turbulence there present. We 
may remark that the values of R 5 are not inconsistent 
with those found in flow in pipes. The speeds ob- 
tained from the pressure distribution are undoubtedly 
of reasonable accuracy, but the assumption of a linear 
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Figure 24.— Drag coefficient for wooden replica of II-Q-12 model at an 
air speed of 80 ft./sec. with wire rings of varying diameter and position 
on the model. x= distance of ring aft of nose 

distribution in the boundary layer gives too small a 
thickness. In the case of skin friction on plates, equa- 
tion (10) shows that this approximation gives a thick- 
ness only 63 per cent of the value obtained by the pre- 
cise computation. Hence, instead of 2,320 as the low 
value of the Reynolds Number, we have a smaller value. 

These curves should not be taken too seriously, al- 
though the agreement is better than could reasonably 
be expected. The assumption of discontinuous changes 
in the type of flow, the omission of the form resistance 
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Figure 25— Comparison of computed and observed drag coefficients on N. P. L. short model 
Curves (J , 4H-foot tunnel, upstream section, R h for transition assumed =1,250. 
Curves K , 4^-foot tunnel, working section, R h for transition assumed =2,000. 
Curves D, 4^-foot tunnel, downstream section, R s for transition assumed =2,750. 

Turbulence due to diminishing cross section assumed introduced at x=lA feet. 
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Figure 26— Comparison of computed and observed drag coefficients on N. P. L. long model 
Curves (J , 4J/Hoot tunnel, upstream section, R d for transition assumed =1,250. 
Curves K , 4 , / 2-foot tunnel, working section, R b for transition assumed =2,000. 
Curves D, 4^-foot tunnel, downstream section, i? 5 for transition assumed =2,750. 

Turbulence due to diminishing cross section assumed introduced at x=1.7 feet. 
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that may be present to some extent, and the crude 
approximations that have been made are not intended 
to be accurate pictures of the actual phenomena, which 
are extremely complicated. We hope that the discus- 
sion has been suggestive and we believe that the major 
features are in some degree correct. 

REMARKS ON AIRFOILS 

The aeronautical engineer will naturally ask why 
no experiments were made on airfoils. The reason is 
threefold. First, experiments in the same wind tun- 
nels showing large effects on airship models showed 
small effects on an airfoil model. Second, from the 
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Figure 27.— Comparison of computed and observed drag coefficients on wooden 

II-Q-12 model 

Curves (J , 4^-foot tunnel, upstream section, R 5 for transition assumed ==1,250 
Curves K t 4^-foot tunnel, working section, 7? 5 for transition assumed = 2,000. 
Curves D t 4H-foot tunnel, downstream section, It s lor transition assumed =2, 
Turbulence due to diminishing cross section assumed introduced at 

discussion on airship models, it will be seen that small 
effects are to be expected in an ordinary wind tunnel. 
While Figure 15 was computed for the N. P. L. long 
model, curves of the same general nature were found 
for other models for which the pressure distribution 
curves differed considerably. For airfoils the Rey- 
nolds Numbers are small because of the short length. 
Furthermore the skin friction is only a part of the total 
drag. Third, the experimental difficulties of dis- 
tinguishing small effects when testing at several loca- 
tions and in wind tunnels of different size with port- 
able balances are very great. It may be expected, 
however, that at larger Reynolds Numbers, for ex- 
ample, those obtained in the variable density wind 
tunnel, the effects of turbulence might be distinguished. 
The general nature of the effect is easily predicted. 
We might expect first to find differences in the mini- 
mum drag coefficient. In the new variable density 
tunnel, which has less turbulence, lower coefficients 



are to be expected with the scale effect curve showing 
a minim u in followed by a region of increasing coeffi- 
cient. Further, since separation is delayed by in- 
creased turbulence, a somewhat smaller maximum lift 
coefficient may be expected in the new variable den- 
sity tunnel. We do not know whether any such 
effects have actually been found. It may prove that 
they are entirely negligible. 

CONCLUSION 

STATUS OF WIND TUNNEL STANDARDIZATION 

It will now be appreciated that wind tunnels can not 
be standardized in the sense originally intended. It 
is not possible to determine one or more correc- 
tion factors by means of which results on a new 
model may be corrected to be comparable with 
the results of some standard tunnel. It is pos- 
sible, we think, to assign a characteristic number 
to each tunnel, such that wind tunnels having 
the same characteristic number will give compar- 
able results. This characteristic number will be 
the measured turbulence or, more conveniently, 
the Reynolds Number for which the sphere drag 
coefficient is 0.3. 

The only real standardization that could be 
made would be accomplished by insisting that 
all wind tunnels be constructed so as to have the 
same turbulence. Strange to say there is some 
difficulty in agreeing on the ideal amount of tur- 
bulence. It appears to us self-evident from the 
scientific point of view that the ideal is zero tur- 
bulence. But the practical engineer replies that 
the curves of Figure 15 for a tunnel of small tur- 
bulence can not be extrapolated from the usual 
model range to give the value at a high Reynolds 
Number, whereas the curves for very turbulent 
wind tunnels can, at least approximately. Ex- 
pressed physically, the flow about the model in a 
turbulent wind tunnel at low Reynolds Numbers 
is more like the flow about the model at high 
Reynolds Numbers in a nonturbulent stream 
than is the flow in a nonturbulent tunnel at low 
Reynolds Numbers. One answer is a variable den- 
sity tunnel of low turbulence. Another, less satis- 
factory, is the judicious use of wire rings on the model 
to stimulate artificial turbulence, a controllable proc- 
ess in the nonturbulent wind tunnel as compared to 
the use of the turbulent wind tunnel, in which the 
turbulence can not readily be reduced. 

We conclude by stating that turbulence is a variable 
of some importance at all times and that the careful 
experimenter will desire to measure and state its value 
in order that his experiments may be capable of 
interpretation. 
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APPENDIX 



MODIFICATIONS OF APPARATUS FOR MEASURING 
TURBULENCE 

Since the publication of Technical Report 320, the 
apparatus there described for the measurement of 
turbulence has undergone several important modifica- 
tions, which we wish to describe. We have called 
attention in Technical Report 320 to the fact that the 
calibration of the hot-wire anemometer was very 
unstable, the calibration often changing so much 
during the course of an afternoon that all observa- 
tions had to be discarded. We now believe that a 
large part of the instability was due to the use of soft 
solder for attaching the wire to its holder. We believe 
that the wire was held mechanically without intimate 
contact so that contact differences of potential were 
always present. At any rate we find a spot welding 
method used by our atomic physics section much more 
satisfactory in that the changes in calibration of the 
wire with time are very markedly reduced. The spot 
welding is done by holding the wire against its steel 
support by a copper electrode and momentarily passing 
the short circuit current of a stepdown transformer 
(110-volt primary, 11 to 1 ratio, 1 Kva rating) through 
the electrode and support. The tapping of a key in 
the primary circuit sends a rush of current through the 
contact between copper and steel which develops a 
temperature great enough to melt a little iron around 
the platinum wire. The copper electrode does not 
melt and stick to the wire because of its greater heat 
conductivity. 

The second modification introduced is the use of a 
12-volt heating battery instead of the 120-volt battery 
line, an absolutely essential modification if the appa- 
ratus is to be at all portable. The effect of this change 
is to make the fluctuation of the heating current during 
the speed fluctuations of appreciable magnitude so 
that the calibration curves for constant heating cur- 
rent can no longer be applied directly. It is neces- 
sary to modify the method of computing the speed 
fluctuation from the observed voltage fluctuation. As 
we are most interested in small fluctuations such as 



occur in wind-tunnel air streams in the absence of a 
model, we may consider the calibration curve linear 
over the interval in question and use the process of 
differentiation. The calibration curve of the wire 
according to equation (3) of Technical Report 320 is 

K — iio 

where R 0 is the resistance of the wire at room tempera- 
ture, R is the resistance of the wire when heated in 
an air stream of speed, C7 0 , by a current, i, a is the tem- 
perature coefficient of resistance, and K and C are 
constants. Differentiating this expression, permitting 
% and R to vary, we have 



2i R R 0 a 



(1/ 



Ro 2 



a dR = CdUo 



R — Rq u * (i?-i? 0 ) 2UXt 2Vc/o 
To connect di and dJ?, we have the relation 
12=; (B + r) 

where 12 is the battery voltage and r is the resistance 
of the heating circuit, excluding that of the wire. 

Hence m _ 

i 9n -i 2 dR 



di- 



R + r 



d/r 



12 



and we find on substitution, setting idR = dPJ, the 
measured voltage fluctuation, an approximation which 
is very close: 

dU 0 _ 2 f iR 0 2 a , lVRRoaljj? 

^""^oU^-^o^e R-Ro _T * 

The second term in the bracket represents the correc- 
tion for the variation of the current. 

A typical run at the working section is given in 

Table IX. C is obtained from the plot of % vs 

Vt/o (not shown) as 0.000151. All computations are 
made by slide rule with sufficient accuracy. It is 
seen that the correction for the current variation is 
from 5 to 15 per cent. 
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TABLE IX.— METHOD OF COMPUTING SPEED FLUCTUATION FROM VOLTAGE FLUCTUATION 

[i?o= 3.750 ohms, a =0.0037, resistance of leads to wire, 0.527 ohms, mean heating current, 0.2 ampere] 



Air 
speed, 

U* 
ft./sec. 


Vt/b 


E 
(volts) 


/?+leads 
ohms 


R 
ohms 


R-Ro 
ohms 


A 


B 


D 




F 


<1E 


dUo 
Uo 


29.0 
30.2 
37.0 
45.4 
52.4 
61.4 
67.2 
72.6 

Mean. 


4.90 
5.50 
6. 08 
6. 73 
7.23 
7.83 
s. 20 
8. 52 


1.432 
1.366 
1.315 
1.267 
1.240 
1.214 
1.195 
1. 186 


7.16 
6.83 
6. 575 
6. 335 
6.20 
6. 07 
5. 975 
5. 93 


6. 633 
6.303 
6. 048 
5.808 
5. 673 
5. 543 
5. 448 
5. 403 


2. 883 
2. 553 
2. 298 
2. 058 
1.923 
1.793 
1.698 
1.653 


0. 0001922 
2172 
2414 
2696 
2888 
3094 
3267 
3256 


0. 001274 
1369 
1460 
1566 
1635 
1715 
1779 
1814 


0. 001250 
1596 
1970 
2453 
2812 
3230 
3607 
3805 


0. 001462 
1824 
2213 
2714 
3084 
3516 
3903 
4107 


0. 395 
.439 
.483 
.534 
.565 
.595 
.631 
.638 


-0.00411 | 

350 
316 
296 

2*6 
253 
243 
233 


0. 0162 
.0154 
.0153 
. Oi:>.K 
.0162 
.0151 
.0153 
.0149 
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The third modification is the use of a resistance- 
con denser-coupled amplifier instead of the direct 
current amplifier. With the direct current amplifier 
considerable difficulty was experienced from drift and 
consequent change of amplification factor owing 
to operation at different points on the tube character- 
istic curves. By using very large condensers as 
coupling condensers, it is possible to pass very low 
frequencies. We have used high grade mica con- 
densers, each condenser having a capacity of 2 micro- 
farads. In conjunction with 1 megohm grid leaks, 
the time constant of the coupling circuit is 2 seconds. 
We shall state without proof that the computed errors 
for a frequency of 1 cycle per second are of the order 
of 1 per cent in amplitude and 5° in phase. 

The operation of this amplifier has been found very 
satisfactory. Adjustments arc required infrequently 
and the measured amplification factor remains constant 
for long periods. It is of course necessary to use 
alternating current for measuring the amplification. 
We have been able to use a single reversing commu- 
tator by replacing the inductance of the compen- 
sating circuit by an equivalent resistance. The 
same alternating current instrument may then be 
used in the potentiometer circuit and in the output 
circuit so that the errors due to the rectangular wave 
shape are small. 

One feature of the amplifier which is occasionally 
troublesome is the fairly long time required for the 
effect of transient disturbances to disappear, i. e. the 
effect of adjustments, of a gust of wind in the outdoor 
tunnel or of a variation in the line voltage. It is 
only on rare occasions that the transients are so fre- 
quent that measurements are not possible. 

We may state again the frequency range covered. 
The lower limit for errors of the order of 1 per cent is 
about 1 cycle per second. The upper limit is fixed by 
the functioning of the compensating circuit and is 



about 100 cycles per second. While even at 500 
cycles per second the errors computed from the theory 
of the compensating circuit are not large, the effect 
of the distributed capacity of the large inductance in 
the compensating circuit enters to increase the error 
and we therefore state the useful frequency range as 
1 to 100 cycles per second. 

The fourth modification of our apparatus is the sim- 
plification of the accessory circuits and of the general 
arrangement. Figure 28 shows the wiring diagram 
of the modified arrangement. A photograph has 
already been given in Figure 3. We have made use of 
jacks and plugs to simplify the switching arrange- 
ments. Beginning at the upper right-hand corner of 
the wiring diagram (fig. 28) we have a standard cell 
connected to an open jack. To the left of this circuit 
is the heating circuit. Potential leads from the wire 
and from the manganin resistance used for accurate 
measurement of the heating current are taken to 
open jacks. To the left of the heating circuit is the 
potentiometer circuit, the balancing circuit of which 
ends in a plug. The potentiometer plug may be placed 
in the standard cell jack to measure the voltage of the 
potentiometer battery; in the hot-wire jack to measure 
the mean voltage drop across the wire; into the current- 
measuring jack for determination of the heating cur- 
rent; or into the input jack of the amplifier to calibrate 
the amplifier. The circuit consisting of two plugs and 
a closed jack at the extreme left is used to pass on the 
voltage fluctuations to the amplifier. After the mean 
voltage drop has been balanced by plugging the poten- 
tiometer into the hot-wire jack, the potentiometer 
plug is withdrawn and placed in the jack of the auxil- 
iary circuit. One plug of the auxiliary circuit is placed 
in the hot-wire jack, the other into the input jack of 
the amplifier. The connections are so arranged thai 
only the fluctuations of the voltage about the mean 
value are impressed on the amplifier. 
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The amplifier circuit is shown below the accessory 
circuits in Figure 28. Attention is called only to the 
special features, namely, the use of separate A and B 
batteries for the power stage, the subdivision of the 
plate resistance in the first stage for varying the ampli- 
fication by a factor of one-fourth, one-half, or three- 
fourths, and the special compensating circuit. All 
of these features are discussed in Technical Report 320. 
There are omitted from the wiring diagram several 
battery switches, jacks in each plate circuit by means 
of which plate currents may be checked and stages 
omitted to reduce the amplification, the plugs and 
jacks between the compensating inductance and the 
amplifier and between the milliammeters in the out- 
put circuit and the amplifier, and a switch for removing 
the inductance during calibration of the amplifier. 
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We believe that with these modifications results are 
more accurate, and we know that operation is more 
convenient. While Table IX shows a variation from 

the mean value of of only 0.0006, we have ob- 
served larger differences and for the present can not 
guarantee values to better than 0.002 as previously 
stated. 
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LIST OF SYMBOLS 

= tangential component of the velocity of the 

fluid at any point in the boundary layer. 
= normal component of the velocity of the fluid 

at any point in the boundary layer. 
= speed of the fluid at the outer edge of the 

boundary layer. 
= speed of the fluid at a great distance, i. e., the 

wind tunnel speed. 
-U— 2^ = speed of the fluid at any point in the 

boundary layer relative to the speed at the 

outer edge. 

: distance measured from the leading edge or 
nose along the surface in a plane parallel to 
the wind direction. 

: distance measured normal to the surface. 

= thickness of boundary layer. 

s radius. 

= diameter. 

-- overall length of body. 

: ^j— = maximum area of cross section taken nor- 
mal to the wind direction, 
total surface area from the nose to a distance 

x from the nose, 
force per unit breadth across the stream, 
force per unit area at any point. 
: average force per unit area, 
static pressure, 
density of the fluid. 
Force 



v- 

li- 
lt, > 
C F = 



: viscosity of the fluid. 
: m/p = Kinematic viscosity 

m 



U5 



= Reynolds Number of the general flow. 
= Reynolds Number of the boundary layer. 



or 



F a 



h P U 0 2 x 
area. 
Force 



77-2 = force coefficient based on 



=^ = force coefficient based on 



c = 



K= 



1/2 p (Vol) 2 * 3 CV 
volume. 

force coefficients for skin friction on flat plates 
of lengths l r and l 2 in turbulent flow. 

force coefficient for skin friction on that part of 
the surface of a plate of length l 2 (the overall 
length of the plate) between x — b and the 
rear edge. 



1 rvJZY 

fuy \Uo) 
\u 0 J Jo 
/; w dx 



dx 



R Q = 
R = 

i = 
a = 
r = 



Appendix 

resistance of wire anemometer at room tempera- 
ture. 

resistance of wire anemometer when heated by a 

current, i y in a stream of speed, U 0 . 
heating current. 

temperature coefficient of resistance, 
resistance of the heating circuit excluding that 

of the wire, 
voltage drop across the wire and leads. 
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Positive directions of axes and angles (forces and moments) are shown by arrows 



Axis 


Force 
(parallel 
to axis) 
symbol 


Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


Designa- 
tion 


Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal ___ 
Lateral. _ 


X 
Y 


X 
Y 


rolling.. -__ 


L 
M 


Y i Z 


roll 




u 


V 


pitching 

yawing 


z >x 


pitch___ 


e 


V 


Q 


Normal. 


Z 


Z 


N 


X > Y 


yaw- ___ 




10 


r 















Absolute coefficients of moment 



qbS 



Cm 



q_cS 



Ln ~~ qfS 



Angle of set of control surface (relative to neu- 
tral position), 5. (Indicate surface by proper 
subscript.) 



4. PROPELLER SYMBOLS 



D, 

Pm 



Diameter. 
Effective pitch. 
Mean geometric pitch. 
Standard pitch. 
Zero thrust. 
Zero torque. 
p/D, Pitch ratio. 
V Inflow velocity. 
V 9 , Slip stream velocity. 



!T, Thrust. 
<2, Torque. 
P, Power. 

(If "coefficients" are introduced all 
units used must be consistent.) 
7), Efficiency =T VI P. 
n, Revolutions per sec, r. p. s. 
A 7 , Revolutions per minute, r. p. m. 



<£, Effective helix angle = tan 



5. NUMERICAL RELATIONS 



1 hp = 76.04 kg/m/s = 550 lb./ft./sec. 
1 kg/m/s = 0.01315 hp 
1 mi./hr. = 0.44704 m/s 
1 m/s = 2.23693 mi./hr. 



1 lb. =0.4535924277 kg 
1 kg = 2.2046224 lb. 
1 mi. = 1609.35 m = 5280 ft. 
1 m = 3.2808333 ft. 



